Reaction mechanisms for two-neutron halo breakup 
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We describe a reaction mechanism which is consistent with all available experimental information 
of high energy three-body breakup processes. The dominating channels are removal of one of the 
three halo particles leaving the other two either undisturbed or absorbed. We compare with the 
commonly used deceptive assumption of a decay through two-body resonance states. Our predictions 
can be tested by measuring neutron-neutron invariant mass spectra. 
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Introduction. Halo states are basically characterized 
as spatially extended weakly bound systems. Two-body 
halos interacting with a target present a three-body scat- 
tering problem, which in practice is much more com- 
plicated due to the intrinsic structure of the three con- 
stituent particles. Three-body halos present analogously 
at least a four-body problem, which has to be approx- 
imated preferentially by using physical insights. Reac- 
tions for high beam energies allow a separation of the 
degrees of freedom related to the fast relative projectile- 
target coordinates and the slow intrinsic halo motion. 

The properties of halo systems have been discussed 
intensely over the last decade both in dripline nuclei 
[Q-Q and in molecular systems (|||. Precise definitions, 
classification and occurrence conditions were recently at- 
tempted although different from an earlier defini- 
tion M. Few-body concepts and techniques are success- 
fully applied in the descriptions JI(J , which to a large ex- 
tent focussed on nuclear three-body halos. Much efforts 
have been devoted to two-neutron Borromean halos like 
6 He (n+n+ 4 He) and n Li (n+n+ 9 Li) where two neutrons 
surround a core pj . The basic structure is essentially 
agreed upon while reaction descriptions and analyses of 
measurements still are controversial. 

Halo physics is a substantial part of experimental pro- 
grams with radioactive beams and it is urgent to root out 
widespread misconceptions and clarify how the reactions 
proceed. Furthermore these questions are of general in- 
terest as basic few-body reaction problems. Since the 
halo concept now is applied and exploited in molecular 
physics we also may anticipate similar implications of 
properly formulated reaction models. 

The purpose of this letter is to (i) establish the re- 
action mechanism for two-neutron halo breakup in high 
energy collisions with light targets and in passing clarify 
the differences to the entirely different mechanism due 
to the large charges of heavy targets, (ii) investigate the 
validity of the erroneous but commonly used assumption 
of breakup through resonances in the two-body subsys- 
tems. These questions are crucial and answers urgently 
needed for understanding reactions with halo nuclei. 



Reaction mechanisms. The dominating reaction 
channels for two-neutron halo breakup on light targets 
are experimentally established (^[lSl and theoretically 
described 0,|4| as removal of one neutron or destruc- 
tion of the core, thereby leaving the final state with the 
core and the other neutron or with the two neutrons. 

The decisive question in this context is which reac- 
tion mechanism is responsible for the observed behav- 
ior? The reaction time for light targets is short compared 
to the time scale of the intrinsic halo motion. For spa- 
tially extended systems the target can then remove one 
of the halo particles instantaneously without disturbing 
the motion of the other two particles. This means that 
the sudden approximation basically is valid as accepted 
in several previous publications [ p^| , ^3|p^Jl^ | . 

The implication is that the remaining two-body system 
is left in its initial state which, as unbound for Borromean 
systems, falls apart influenced by the corresponding two- 
body interaction. This decaying two-body system is thus 
formed as a wave packet consisting of those parts of the 
relative two-body wave function present within the origi- 
nal three-body system, which precisely lead to the domi- 
nating reaction products |f0(| . The surviving wave packet 
then has a large component describing the tail of the two- 
body wave function. The short distance parts lead to a 
large extent to removal of more than one particle at a 
time. All other breakup reactions are analogously de- 
scribed in this participant-spectator model (PSM). 

R-matrix formulation. The observed invariant two- 
body mass spectra and the momentum distributions are 
routinely analyzed as arising from the decays of low-lying 
two-body resonances or virtual s-states |r^ , |f3| , |i5| p~7|| . 
These assumptions are in direct contradiction to the 
short reaction time and the sudden approximation. 
There is not sufficient time for the remaining two halo 
particles to adjust their relative motion and populate 
corresponding resonance states. This requires at least a 
reaction time comparable to the intrinsic halo time scale. 

Thus these analyses apparently invoke both the sudden 
approximation and decay through resonances or virtual 
s-states. These assumptions are strictly incompatible ex- 
cept when these two-body states are populated within the 
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initial three-body system. This is clearly seen by con- 
structing a Borromean system by adding a neutron to 
a neutron-core resonance state. The overlap of this and 
the real bound state wave function may still be substan- 
tial, but rearrangements are necessary to reach the bound 
three-body state, i.e. a novel few-body system carrying 
otherwise inaccessible information about the off-shell be- 
havior of the nucleon-nucleon interaction. 

The pertinent questions are what we can learn from 
measured two-body invariant spectra and what infor- 
mation is in fact obtained by the analyses using reso- 
nances or virtual s-states as intermediate states. The 
analyses, often erroneously claiming to use Breit-Wigner 
shapes E|Jl|,||,0, 

are in fact based on i?-matrix the- 
ory 1^ , ^91 , where a complete basis of two-body contin- 
uum states are used after removal of one of the three 
particles. This basis could consist of "correct" low-lying 
resonance states supplemented with a discretized or con- 
tinuum higher lying set of states. 

In practice one proceeds by reducing the unspecified 
and unknown basis to a few terms, i.e. usually one or two 
states. This reduction of model space may be allowed if 
the basis consistently is renormalized, i.e. the new basis 
includes properties of the excluded states. Thus fitting 
in this context by use of a small basis seems to prohibit 
interpretation in terms of the "correct" two-body reso- 
nance states. In principle maintaining the basis without 
renormalization would be correct but this presupposes 
exactly that knowledge about these states, which is the 
very aim of the analyses. This problem cannot be solved 
by increasing the employed model space until conver- 
gence is reached and no renormalization is needed. A 
larger space implies more parameters in the fitting pro- 
cedure and reproduction of the data is not unique. The 
problem becomes overdetermined and the extracted pa- 
rameters inaccurate or directly unreliable. 

The analyses using two-body resonances or virtual s- 
states therefore assume (i) that no renormalization due 
to truncation of model space is needed, (ii) a reac- 
tion mechanism where only the "clean" two-body res- 
onance states or virtual s-states are populated and (iii) 
no other (known or unexpected) reaction channel con- 
tribute. These assumptions are at least inaccurate. The 
difficulties are enlarged when more than one resonance or 
more than one reaction channel contribute. When the as- 
sumptions are fairly well fulfilled the interpretation would 
also be approximately correct. 

Computations. We shall concentrate on 6 He and n Li 
colliding with light targets. We shall use the PSM formu- 
lation, where one halo particle (the participant) interacts 
with the target while the other two halo particles (the 
spectators) are left undisturbed [jn^. The participant- 
target interaction is described by the phenomenologi- 
cal optical model while the spectators are treated in 
the black sphere approximation, i.e. they are absorbed 
within a given radius from the target and otherwise they 
continue undisturbed. This model faithfully exploits the 
consequences of a short reaction time. We compute the 



population of the two-body continuum states after the 
instantaneous removal of the third particle. 

The i?-matrix expressions of the invariant mass spec- 
trum da/dE of the spectator system and the relative 
spectator momentum distribution Pi ong are ]T^Jl5| , |l7t 

do _ gj T(E) £'+°- 5 

dE 2ir (E-E r ) 2 +0.25r 2 (£)' £<+°- 5 ' {> 

f°° dE da 2 . . . 

PlongiP) = J ^/EdE 7 E min=P / 2 M ( 2 ) 

where <j\ is the total cross section, I is the orbital an- 
gular momentum, /i is the reduced mass, E r and Tq are 
position and width parameters. The distributions are 
correlated and should not be fitted independently. Pre- 
cisely the same procedure applies both when the core and 
a neutron are the participant, i.e. the final state consists 
of two neutrons or a neutron-core system, respectively. 

The chosen observables amplify the effects of the as- 
sumed reaction mechanism. We can then compare the 
experimental distributions both with the PSM predic- 
tions jlO| and the i?-matrix results obtained by the decay 
through resonance assumption. This provides evidence 
about the basic reaction mechanism. 

Neutron removal. Absorption of one neutron from 
6 He produces a neutron- 4 He continuum state, which 
mainly is of p 3 / 2 -character, since the P1/2 neutron-core 
state has a higher energy and the Si/2-wave is repul- 
sive. The corresponding invariant mass spectrum and 
the relative momentum distribution are shown in Fig. [I]. 
The PSM computation agrees fairly well with the mea- 
sured invariant mass spectrum The peak position 
reflects the energy of 0.77 MeV of the neutron- 4 He P3/2- 
resonance with the width of 0.5 MeV used in the PSM 
computation. Any value of r , see Eq.(|), from 0.4 MeV 
to 0.8 MeV also reproduce the experiment fairly well. 

For consistency the momentum distribution should 
now follow with the same parameters. Indeed we see in 
Fig. [l] that the PSM results resemble the (almost) width 
independent i?-matrix fits confirming that the width pa- 
rameter for the n- 4 He resonance can not be determined 
in this way. We also note the characteristic flat maximum 
of p— waves. 

The 11 Li system is different due to the core spin of 
3/2 and the mixture of s and p- waves in the subsys- 
tems. The computed three-body wave function contains 
around 60% of s-wave and 40% of p 2 -wave neutron-core 
configurations. The neutron- 9 Li system has a low lying 
virtual s-state at 240 keV and a p-resonance at 0.5 MeV 
ifiof . Neutron removal results in the distributions shown 
in Fig. |^. The contribution to the invariant mass spec- 
trum from s-waves peaks at a very low energy determined 
entirely from the phase space constraint and independent 
of the position of the virtual s-state jl9| . In contrast the 
p-wave contribution peaks at the two-body resonance en- 
ergy. The measured spectrum is fairly well reproduced 
by the PSM computation again supporting the assumed 
initial three-body structure and the reaction mechanism. 
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FIG. 1. Neutron- 4 He invariant mass spectrum (upper) 
and longitudinal relative momentum distribution (lower) for 
breakup of 300 MeV/nucleon 6 He projectile on a carbon tar- 
get. Experimental data from pl| . The solid curves are the 
PSM calculations and the broken lines are obtained from 
Eqs.(^) and ([2J) with 1 = 1. The invariant mass curves have 
been convoluted with the instrumental response po| . 

We compare in Fig. || with two different i?-matrix fits. 
In the first the computed s and p contributions are fitted 
separately thereby maintaining the samep-wave content. 
In the second fit we use the parameters in , which also 
reproduces rather well the experimental (and PSM) data. 
The 10 Li structures underlying these fits differ substan- 
tially as expressed clearly through the different widths. 
The p-wave contents also differ substantially, i.e. about 
35% for the first and 70% for the second fit @. 

In the lower part of Fig. |^ we show the corresponding 
relative momentum distribution. The PSM computation 
and the first fit produce a very similar momentum distri- 
bution, while the second fit differs in the central part due 
to the large fraction of p-waves that creates the plateau 
at low relative momentum. Therefore different fits of 
invariant mass spectra of similar accuracy can produce 
rather different momentum distributions due to empha- 
sis of different features of the distribution. 

Core destruction. The reaction assumptions can be 
tested by similar investigations of the other important 
channel corresponding to destruction of the core and 
leaving the two neutrons as spectators. The final states 
then consists of two neutrons for both 6 He and 11 Li. Thus 
we can separate effects of initial and final state structures. 
However, this assumes that fragments from the core in- 
teracting strongly with neutrons are excluded from the 



data. In the PSM computations, discussed in connection 
with Figs. [l] and 0, the spectra are sensitive to the initial 
three-body structure. The root mean square distance 
between the neutrons is more than 6 fm for 11 Li and less 
than 4.5 fm for 6 He. The neutron-neutron invariant mass 
spectrum and the corresponding momentum distribution 
are then both expected to be substantially narrower for 
11 Li than for 6 He. The same consistent PSM model has 
been tested on many other, relative and absolute val- 
ues of observables for neutron removal and core breakup 
reactions for both projectiles |l(|. The agreement with 
available experimental data is overall very convincing. 
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FIG. 2. The same as Fig. [j] for neutron- 9 Li for a "Li pro- 
jectile. Both sets of experimental data p^|p7| are normalized 
to the same maximum as the computed spectrum. The thick 
solid curves are the PSM calculations consisting of contribu- 
tions from both neutron-core relative s (thick dashed) and p 
(dotted) waves. The thin solid and dashed lines are weighted 
averages of two distributions from Eqs.(jjj) and (^) correspond- 
ing to s and p-resonances with parameters given in MeV. 

Further tests of the PSM model (and the i?-matrix 
analyses) would be measurements comparing to the pre- 
dictions presented in Fig. ^. The neutron-neutron rela- 
tive s- waves are completely dominating for both 6 He and 
11 Li and consequently the invariant mass spectra have 
very low-lying peaks. Both spectra and momentum dis- 
tributions are qualitatively similar for the two cases, but 
quantitatively the 11 Li results are much narrower than 
those of 6 He. The i?-matrix distributions fitting the two 
PSM curves in Fig. || correspond to very different en- 
ergy and width parameters without any connection to the 
known neutron-neutron scattering properties. The PSM 
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model predicts different neutron-neutron spectra for 6 He 
and n Li after core breakup. A reaction mechanism popu- 
lating final state two-body resonances independent of the 
initial structure must predict identical neutron-neutron 
invariant mass spectra for both projectiles. Experimental 
data could distinguish between these models. 
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FIG. 3. The same as Fig. [ljfor the neutron-neutron system 
for 6 He and lx Li projectiles. The points are the two _R-matrix 
fits given in the figure. 

Conclusions. The dominating reaction channels for 
high energy breakup of Borromean three-body halos on 
light targets are one-particle removal and subsequent de- 
cays of the wave packets created in these processes. The 
reaction time is short and any resonance structure of 
the remaining two-body system is populated with the 
amount already present in the initial three-body wave 
function. All available experimental data for high energy 
three-body breakup on light nuclei are consistent with 
this reaction mechanism. For heavy targets the reaction 
mechanism for the dominating channel is quite different 
proceeding through a gentle excitation of the three-body 
continuum by the Coulomb interaction. 

Analyses assuming instantaneous removal of either a 
neutron or the core, while populating resonances in the 
remaining two-body system, are conceptually inconsis- 
tent. An invariant mass spectrum reproducing the data 
only reflects that the corresponding energy distribution 
was present immediately after the final state two-body 
system was isolated. The problem is especially enlarged 
when more than one resonance or virtual state are impor- 
tant for the two-body subsystems. The inconsistency is 
highlighted in spectra obtained after core breakup, where 



the final states are identical (two neutrons). Then the 
resulting distributions should also be identical for dif- 
ferent two-neutron halo projectiles, even when the initial 
three-body structure differs substantially. This is in clear 
disagreement with elaborated consistent model compu- 
tations reproducing essentially all available data. In any 
case, the neutron-neutron invariant mass spectra provide 
direct evidence of the breakup reaction mechanism. 

We conclude that correct interpretation of the invari- 
ant spectra almost inevitably require a consistent model, 
i.e. the three-body structure, the two-body interactions 
and the reaction mechanisms must be (approximately) 
correct. Even then interplay between the different in- 
gredients may produce misleading results. The i?-matrix 
analyses rely on assumptions or computations of the ini- 
tial steps producing an isolated two-body system. 
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